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Abstract: The diffraction trace formula (Phys. Rev. Lett. 73, 2304 (1994)) 
and spectral determinant are tested on the open three disk scattering sys- 
tem. The system contains a generic and exponentially growing number of 
diffraction periodic orbits. In spite of this it is shown that even the scatter- 
ing resonances with large imaginary part can be reproduced semiclassicaly. 
The non-trivial interplay of the diffraction periodic orbits with the usual ge- 
ometrical orbits produces the fine structure of the complicated spectrum of 
scattering resonances, which are beyond the resolution of the conventional 
periodic orbit theory. 
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1 Introduction 



Gutzwiller's trace formulao is an increasingly popular tool for analyzing 
semiclassical behaviour. Recently, it has been demonstrated that using proper 
mathematical apparatus, like the spectral determinant of Vorosa, cycle 
expansions® or quantum Fredholm determinants!, the trace formula can 
successfully predict individual eigenenergies of bound systems and resonances 
of open scattering systems. The physical content of the trace formula is the 
geometrical optical approximation of quantum mechanics via canonical in- 
variants of closed classical orbits. This approximation is very accurate when 
periodic orbits sufficiently cover the phase space of the chaotic system. This 
is not the case when the number of obstacles is small or their distance is 
large compared to their typical size. Such a problem occurs where the wave 
length of a quantum mechanical (or optical) wave is very large compared 
to the spatial variation of a repulsive potential, e.g. at the boundaries of 
microwave guides, optical fibers, superconducting squids, or circuits in the 
ballistic electron transport, i.e. in most of the devices used for so-called meso- 
scopic physics. In such cases it is important to take into account the next- 
to-geometrical effects. In Ref.|| we have shown how the Geometric Theory 
of Diffraction (GTD) for hard core potentials can be incorporated in the 
periodic orbit theory. We worked out the two disk scattering system as an 
example, where diffraction plays an important role. Since the realization of 
the importance of such effects, diffraction periods have been uncovered in 
rhomboid billiardsB , billiards with magnetic flux lines@ and in a limiting 
case of the hyperbola billiardcS In the present example we study for the first 
time a generic example where an exponential number of diffractive and geo- 
metrical orbits interplay and builds up a complicated spectrum of scattering 
resonances. For the reader unfamiliar with the diffraction trace formula we 
start by a brief sketch of its derivation. For a more detailled and completely 
outlined introduction to the theory we refer the reader to [?]. 



2 Diffraction periodic orbits 

As is well known the free particle Green function G(x,y; E) can be exactly 
described in terms of geometrical optics by the path that connects x with y 
at energy E. If a smooth potential is introduced we have to respond with 
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a refractive index and if hard walls are present have to deal with diffractive 
rays to keep the description in the spirit of geometrical optics. The diffractive 
rays connecting two points in the configuration space can be derived from 
an extension of Fermat's variational principle of classical mechanics : Each 
path connecting x with y has a whole class of topologically equivalent paths 
which can all be continuously deformed into each other without changing 
the number of encounters with the hard wall singularities of the system. The 
generalized Fermat principle then states that for each such class T only the 
rays of stationary optical length among all the curves in T contributes to the 
final field. The total field will then be the sum of contributions from such 
paths with diffractive segments, over all the topologicaly different classes of 
paths connecting x with y. 

Once we know the generalized ray connecting two points A and B we can 
compute semiclassicaly the Green function (fc, E) by tracing the ray@: 

a, On the geometrical segments of the ray, the Green function is given by 
the energy domain Van-Vleck propagator 

G(q, Q ', E) = j^r 2 D^( q , q', E)e^'^, (1) 

where Dy(q,q', E) = | det(— d 2 S/dqidq'j)\/\q\\q'\ is the Van-Vleck determi- 
nant and v is the Maslov index (see Ref. [|1(| for details). 

b, When the geometrical ray hits a surface, an edge or a vertex of the 
obstacle it creates a source for the diffracted wave. The strength of the source 
is proportional to the Green function at the incidence of the ray 

Q diff = .DGinci • (2) 

The diffraction constant D depends on the local geometry of the obstacle, 
the wave-type and the nature of the diffraction. It has been determined in 
Ref. |9j from the asymptotic semiclassical expansion of the exact solution 
in a simple geometryB Q. For the surface diffraction (creeping) it has the 
form 

D = 2 l/3 3 -2/ 3 W12 (M^! . ( 3 ) 

Here Ai'(x) is the derivative of the Airy integral Ai(x) = dt cos(xt — t 3 ), 
k = V 2mE I Ti is the wave number, p is the radius of the obstacle at the 
source of the creeping ray and xi are the zeroes of the Airy integral. The 
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index I > 1 refers to the possibility of initiating creeping rays with different 
modes, each with its own profile. In practice only the low modes contribute 
to the Green function. For wedge diffraction the diffraction constant is 

I I 

cos(7r/n) — cos((# — ot)/n)) cos(7r /n) — cos((6 l + a + 7r)/n) 

where (2 — n)ir is the angle of the wedge (n is a real number), a is the incident 
angle and 9 is the outgoing angle. For details we refer to Ref. ||. In the 
three disk scatterer the edge diffraction is only important when the system 
is closed and we shall therefore not go deeper into the subject here. 

The source, created by the incident ray, then initiates a new ray propa- 
gating along the surface (for creeping) or a free ray starting at the edge of 
the obstacle (wedge diffraction). 

During the creeping of the ray the amplitude decreases, which can be 
understood as a process analogous to the radiation processes of electrody- 
namics. The radiated intensity is proportional to the intensity of the ray 

^-A l (s,E) 2 = -2a l (s,E)A l (s,E) 2 , (4) 
as 

where s is the length measured along the surface and Ai(s, E) is the complex 
amplitude of the Green's function along the surface. The coefficient ai(s, E) 
depends on the local curvature of the surface, l/p(s), and it has the structure 
oti(s,E) = xie~ l7T / 6 (k/6p(s) 2 y/ 3 (see Ref. fl^l ) , where the index I refers again 
to the different modes of the creeping wave. The Green's function for the 
creeping ray of mode I is then given by 

Gf{q A >,q&,E) = e -^ dsatis ' E) e^ s ^'^'' E \ (5) 

where L is the length of the arc of the creeping ray, and S(q A >, q&, E) is the 
action along it. 

When the creeping ray leaves the surface its intensity can be calculated 
from the relation (^) due to the reversibility of the Green function. The total 
Green function is then the product of the Green functions and diffraction 
coefficients along the ray. If for example we have geometrical propagation 
from A to A! followed by surface creeping from A! to B' and then again a 
geometrical propagation from B' to B, the total semiclassical Green function 
is 

G(QA, Qb) = G(q A , qA>)D A ,G Cr (q A ,, q B >)D B 'G(q B ,, q B ). (6) 



D 



sin(7r/n) 



n 
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Contrary to the pure geometrical case the semiclassical energy-domain 
Green function for rays with diffraction arcs thus have a multiplicative com- 
position law. 

When we incorporate diffraction effects into the trace formula, periodic 
rays with diffraction segments also contribute. We can handle separately the 
pure geometric cycles and the cycles with at least one diffraction arc or edge: 

Tr G(E) « Tr Gq(E) + Tr Gd(E), (7) 

where Tr Gg{E) is the ordinary Gutzwiller trace formula, while Tr Gd(E) is 
the new trace formula corresponding to the non-trivial cycles of the GTD. 
The Gutzwiller trace formula for two dimensional billiards is 

i oo „irS p (E)— irVpir/2 

Tr G G (E) = -VV T P (E)— ? (8) 



where T P (E) is the time, S P (E) is the classical action, v v is the Maslov in- 
dex and A p is the stability eigenvalue of the primitive periodic orbit. The 
summation goes for all primitive periodic orbits of the system p and their 
repetitions r. TtGd(E) can be computed by using appropriate Watson con- 
tour integrals!!!! . For technical details we refer the reader to Refs. [|ll|, |l^ . 
If we denote by i = 1, . . . , n p (with q np+ i = <&) the points along the closed 
cycle, where the ray changes from diffraction to pure geometric evolution 
or vice versa, or where the ray encounters a wedge diffraction, the trace for 
cycles with at least one diffraction arc can be expressed as the product 

1 OO n p 

Tr Gd(E) = tEE T p( E ) Jl[D( qi )G(q u ft+i, E )Y (9) 

m p r =l t=l 

where T P (E) is the time period of the primitive cycle and D(qi) is the diffrac- 
tion constant (JD at the point qi. G(qi, qi+i, E) is either the Van-Vleck prop- 
agator, if qi and g^+i are connected by pure geometric arcs, or is given by 
the creeping propagator (|5]) in case q^ and q i+ i are the boundary points of a 
creeping arc. 



3 The 3-disk system 

To investigate the theory sketched above we apply it to the 3 disk scattering 
system. The 3-disk system has in recent years been subject to a large number 
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of investigations, and its main virtues are wellknown. Here, we just would 
like to remind about some of the basic properties of the system. 

The system consists of three identical disks placed symmetrically around 
the origin in a plane (see fig. 0), and is completely determined by a single 
parameter, namely the ratio R : a of the seperation of the disks to their 
radius. In our calculations we have kept this ratio fixed at R : a — 6. 

The system thus possesses a 3-fold rotational symmetry around the ori- 
gin and has 3 reflection symmetries around the symmetry lines through the 
center. The system is therefore invariant under the point group Cz v ; and 
instead of considering the system as a whole we can restrict ourselves to the 
fundamental domain^-3. The fundamental domain exactly covers the whole 
system when the elements of the point group are applyed to it. The 3-disk 
system and a version of the fundamental domain are shown in figure [I]. For 
sufficiently large spacing of the disks the system has a complete binary 
symbolic dynamics. All the periodic orbits can be described in terms of the 
alphabet {0, 1} where '0' corresponds to a bounce under which the particle 
returns to its starting disk and '1' corresponds to the bounces where the 
particle continues to the next disk. In the fundamental domain there are 
therefore two fixpoints '1' and '0' corresponding to a triangular, and back- 
and-forward bouncing orbit in the full space. All the geometrical orbits can 
be found via a minimization of the path lengths. If one needs a periodic orbit 
following a definite sequence of n disk bounces, one just has to determine the 
length as a function of the n bouncing positions and then to minimize this 
length. That this indeed gives the right periodic orbit follows from geomet- 
rical optics and Fermats principle: when the light (the particle) follows the 
shortest path (of a given symbolic sequence), it will at the same time obey 
the reflection law. 

The surface of the disk in the fundamental domain can be used as a 
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Poincare surface of section. Establishing the bouncing map as in Ref. JT7 
can thus calculate the stabilities of the cycles. Following the outlined scheme 
we arrive at the results displayed in Table |I[ 

A more detailled description of the 3-disk system and the methods de- 



scribed in this section can be found in e.g. Ref. [17[ . 

In order to apply the GTD to the calculation of semiclassical resonances, 
we also have to account for the diffraction (creeping) orbits of the system. 
To give an overview of the work to be done, we start by counting the number 
of periodic creeping orbits to be evaluated. Because of the symmetry of the 
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system we can assume that the creeping orbit always starts tangentially from 
the (half-) disk in the fundamental domain which we label disk number 1. 
Considering first an orbit with no geometrical bounces we see that it has two 
different disks to go to, and for each each disk two different sides to creep in. 
This makes a total of four diffraction orbits of topological length 1. When 
these are folded back into the fundamental domain we see that two of them 
are self retracing. The two other orbits are tracing the same orbit, but in 
opposite directions. If we consider paths of the particle with m bounces, 
we see that there will be 2 n+1 = 2 m+2 periodic creeping orbits of topological 
order n, as for each one of the m bounces the particle can choose between two 
disks. Thus the number of periodic creeping orbits grows exponentially fast 
with the topological length, n, of the orbit. It is quite astonishing however, 
as we will see later, how few of these orbits are in fact needed to get a 
good description of the scattering resonances (including the ones with large 
imaginary parts). The creeping orbits can be described completely by their 
itinerary \ol\oli . . . a n where the a^s are taken from the alphabet {1,2,3} 
and where we do not allow the repeats . . . 11 . . . , . . . 22 . . . and . . . 33 . . . . 
This description contains a double degeneracy due to the fact that the orbit 
has the choice to creep around the final disk clockwise or anti-clockwise. For 
instance, '123' can represent two different orbits which start from disk 1 in 
the fundamental domain, then hit disk number 2 and finally creep around 
the final disk (3) clockwise or anti-clockwise. 

The restriction that the creeping periodic orbits should start and end 
tangentially on one of the disks simplifies the search procedure for them 
considerably: whereas in the case of geometrical n-bounce cycles one had 
to minimize a function of n bouncing parameters, we here only have one 
parameter in play, namely the angle where the creeping orbit leaves the 
initial disk. Suppose now that we want a specific creeping orbit described by 
a series of disk bounces plus the specification of the final creeping domain 
as above. We then scan through all the angles that leave the first disk in 
the fundamental domain. This gives us an interval of angles where the first 
wanted disk is being hit. We then scan this interval for bounces on the next 
disk in the itinerary and so on. Finally we scan the last obtained interval to 
find the angle under which the ray creeps into the wanted side of the final 
disk. 



8 



4 Cycle expansion of the spectral determi- 
nant 



Having established the data material for the 3-disk system as described above, 
we now report on the more technical part of the actual calculation. 

The resonances can be recovered from the Gutzwiller-Voros spectral determinant! 
A(E), which is related to the trace formula as 

TrG(£) = ^lnA(£). (10) 

The full semiclassical determinant can be written as the formal product of two 
spectral determinants, one corresponding to the pure geometrical, and one to 
the new diffractional cycles: A(E) = A G (E)A D (E), due to the additivity of 
the traces. The product is only formal, since the eigenenergies are not given 
by the zeros of Aq(E) or Ad(E) individually, but have to be calculated from 
a curvature expansion of the combined determinant A(E) itself. 
The geometrical part of the spectral determinant is given by 
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1 p irS p (E)-iri/pir/2 \ 

A G (E) = exp I - V V -— 7 v , (11) 

where the summations are over closed primitive (non-repeating) cycles p and 
their repetitions r. The diffraction part of the spectral determinant is 

/ oo i n p \ 

A D (E) = exp - Y, £ - Ii[D(qf)G(qf, qf +1 , E)\ r , (12) 

\ p r=l r i=l ) 

where the summations are over closed primitive (non-repeating) cycles p and 
their repetitions r. The product of Green functions should be evaluated for 
qf belonging to the primitive cycle p. After summation over r, the spectral 
determinant can be written as 

A D (E) = H(l-t p ) (13) 

p 

with 

rip 

t P = I[D(q!)G(<£,<£ +1 ,E), (14) 

i=l 
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where qf belongs to the primitive cycle p. Here the mode numbers / of 
the diffraction constants and the corresponding summations have been sur- 
pressed for notational simplicity; they can easily be restored as e.g. in the 
final expression (pip . 

We can conclude that the diffractional part Ae>(E) of the spectral de- 
terminant shares some nice features of the periodic orbit expansion of the 



dynamical zeta functions&3, and it can be expanded as 



A D (E) = + • ( 15 ) 

p P ,p' 

The weight ([14]) has the following property which helps in radically reducing 
the number of relevant contributions to the expansion: If two different cycles 
p and p' have at least one common piece in their diffraction arcs, then the two 
cycles can be composed to one longer cycle p+p' and the weight corresponding 
to this longer cycle is the product of the weights of the short cycles 

tp+p> — t p ■ tpt. (16) 

As a consequence, the product of primitive cycles, which have at least one 
common piece in their diffraction arcs, can be reduced in such a way that 
the composite cycles are exactly cancelled in the curvature expansion 

na-y =i-e*6, (ir) 
p b 

where tb are basic primitive orbits which can not be composed from shorter 
primitive orbits. In the case of the desymmetrized 3-disk scatterer this applies 
to all the orbits, and we thus get a zeta function exactly of the form (|17D , 
where the sum is over all the prime periodic creeping orbits. 

To get the free flight part of Ad(E) we first consider the semiclassical 
Green function in free space. This is asymptotically (kR 3> 1) given as 

G »(^' £ »-4(^) 1/2 ^' } ' < 18 > 

where R = \q — q'\. If the ray connecting q and q 1 is reflected once or more 
from the curved hard walls before hitting tangentially one of the surfaces, 
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we can keep track of the change in the amplitude by the help of the Sinai- 
Bunimovich curvatures. For a free flight the Sinai-Bunimovich curvature is 
just the inverse of the travelled distance 

k=-. (19) 

r 

When a hard wall is encountered the curvature changes discontinuously as 

«+ = k_ + — — — — (20) 

COS0 

where k± are the Sinai-Bunimovich curvatures right after and before the 
bounce against the wall and c is the curvature of the reflecting surface at the 
point of incidence, whereas <fi is the angle of incidence. 

By computing the curvatures K{ right after the reflections, and knowing 
the distances U between the z-th and the (i + l)-th points of reflections, 



the factor R in the Green function (18) has to be changed to the effective 



rr 

radius R = l lT£Li(l + ^ K i) where l is the distance between q and the first 
point of reflection along the ray as measured from q, and m is the number 
of reflections from the hard potential walls. The effective radius Rf^, the 
length of the geometrical arc and the length of the diffraction part of 
the first twelwe orbits with creeping sections are listed in Table [D]. To each 
cycle in the list, there is a whole sequence of cycles which wind around the 
disk tow times. For these orbits one has to add 27ramw to the diffraction 
length L®. The diffraction part of the spectral determinant is finally given 
by 



A D (k) = l-]T(-inq 

b,l 

1 



a l/3 e i7r/12 e ifc(Lf +L o ) _ QiL p 



fci/eJ^eff 



(21) 



where Ci = 7r 3 / 2 3 _4 / 3 2 -5 / 6 / Ai'(xi) 2 , ai is the creeping exponent and nif, is 
the number of reflections of orbit b from the disk in the fundamental domain. 
The summation for the windings mw gives the factor 1/(1 — e 2n ( ik - a i) a y 
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5 Numerical results and conclusions 



To evaluate the results of the diffraction extended Gutzwiller-Voros spectral 
determinant, we compare the resonances determined by this, to the reso- 
nances determined just from geometrical orbits and to the exact quantum 
resonances. 

The data are displayed in figures ^ and |3|. As one can see the Gutzwiller 
Voros determinant accounts reasonably well for the leading order of reso- 
nances, whereas it fails for the next series. In figure [5], however, we can 
see that - when a few periodic creeping orbits are introduced - the results 
are qualitatively different, and represent much better the trend of the exact 
quantum resonance data. For instance, one can make a one-to-one identifi- 
cation of the quantum and semiclassical resonances, which is not possible in 
the purely geometrical theory, since in that approximation even the number 
of resonances is wrong. 

The series of subleading resonances also approximately defines the lower 
boundary of the region in which the diffractional spectral determinant still 
has a high accuracy and good convergence properties. This can also be seen 
from fig. m since for small Re k and large Im k we have a relatively larger 
deviation between the exact and creeping resonances. 

The errors of the resonances are originated in two sources. 1. The de- 
scription is semiclassical and therefore we use the Van-Vleck propagator in 
(18) instead of the exact propagator, and the semiclassical approximation of 
the creeping propagator in (5). Also, only the / = creeping modes are used. 
2. Only a restricted number of usual and creeping periodic orbits is available 
instead of infinitely many. 

As mentioned earlier the number of creeping periodic orbits in this system 
increases exponentially with the topological length of the cycles. It would be 
natural to expect that this might destroy the simplicity of the semiclassical 
description. We conclude that this seems not to be the case. As we have 
demonstrated, one only need the basic representatives of the creeping families 
to change the picture of the scattering resonances drastically, in the direction 
of the exact quantum resonances. 
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Figure 1: The full 3-disk system with a copy of the fundamental domain. 
Representatives of the creeping orbits of topological length 1 are displayed 
in full space as well as in the fundamental. 



Figure 2: The exact quantum mechanical resonances (diamonds) and the 
pure geometrical Gutzwiller Voros resonances (crosses) in units of 1/a in 
the complex k plane. The resonances belong to the one-dimensional Ai 
representation of the 3-disk system with R : a — 6 . In the semiclassical 
calculation cycles up to topological length 4 has been used. The leading 
resonances close to the real axis are exactly described by the Gutzwiller Voros 
resonances whereas the subleading semiclassical resonances clearly deviates 
from the exact quantum resonances. 



Figure 3: The exact quantum mechanical (diamonds) and the semiclassical 
(crosses) A\ resonances of the R : a = 6 3-disk system. The resonances 
are calculated by including diffractional creeping orbits up to order 4 in the 
GTD. As in the two disk case an improvement of the approximation is clearly 
visible, especially for the second row of the leading resonances as well as for 
the subleading diffractional ones. In the latter case the qualitative trend 
is clearly reproduced. As discussed above the accuracy of the semiclassical 
resonances becomes worse in the region where Re k is small and Im k is large. 
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Table I: Geometrical cycle data for the 3-disk system with R : a — 6 . The 
first column indicates the symbolic dynamics of the periodic orbit, whereas 
the second and third column gives the stability calculated from the Jacobian 
of the bouncing map, and the length of the cycle in the fundamental domain. 
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Table II: Creeping cycle data for the 3-disk system with R : a — 6. The 
first column indicates the itinerary of the orbit, second column the effective 
radius of the orbit calculated by means of the Sinai-Bunimovich curvatures 
and finaly the third and fourth columns shows the length of the free flight 
and the creeping sections respectively. 
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